Polarized resonant inelastic x-ray scattering as an ultra-fine probe of excited states in 

La 2 Cu0 4 
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X-ray absorption is the standard method to probe the unoccupied density of states at a given 
edge. Here we show that polarized Resonant Inelastic X-Ray Scattering in La2Cu04 at the Cu 
K-edge is extremely sensitive to the environment of the Cu atom and the fine structure in the Cu 4p 
density of states. Combined ab initio and many-body cluster calculations, used for the first time in 
such a context, show remarkable agreement with experiment. In particular we identify a non-local 
effect namely a transition to off-site Cu 3d states. 

PACS numbers: 71.15.Mb, 74.72.Dn, 78.70.Ck, 78.70.Dm 



X-ray absorption spectroscopy (XAS) has long been 
used to investigate the nature of the partial density of 
excited states at a given absorption threshold. However 
pre-edge features, mainly composed of quadrupolar ex- 
citations, are often hard to detect in XAS due to their 
small intensity and the broadening due to the core-hole 
lifetime. As a consequence important insights related to 
the local environment of the absorbing atom are not ac- 
cessible in XAS. Here we show that resonant inelastic 
X-ray scattering (RIXS) yields a precise determination 
of the low-energy unoccupied electronic states and the 
chemical environment of the absorbing atom. RIXS is a 
powerful tool for the study of the electronic excitations 
m solids 1] giving access to site, element and orbital se- 
lective information. It is a second order process where 
the X-ray emission is measured as the incident X-ray en- 
ergy is tuned through an absorption edge of the sample. 
The incident energy determines the intermediate state 
and the core-hole enhances the cross section 

In this work we study resonant Cu ls-2p emission 
in La2CuC>4 both from an experimental and theoretical 
point of view. The excited states in La2Cu04 are of in- 
terest because of the correlated nature of this compound 
which on hole doping becomes a high temperature super- 
conductor. The hopping between the Cu 3d x 2 _ y i and the 
O 2p orbitals combined with Coulomb repulsion between 
Cu 3d electrons form the basis for describing the low en- 
ergy excitations in this solid 0] . In the intermediate state 
of the resonant process a Is Cu core-electron is excited to 
states with energies around the pre-edge region of the Cu 
K-absorption edge. The final state involves Kev emission 
with filling of the Is hole at the expense of a hole in the 2p 
state. Identifying the measured features and explaining 
their dependence on incident energy and incident pho- 
ton polarisation is not a trivial task from the theoretical 



point of view, particularly in strongly correlated systems 
where single particle theories like standard density func- 
tional theory (DFT) dramatically fail due to the lack 
of important many-body effects. A different approach 
to the problem is the use of many-body cluster models 
based on multiband Hubbard-like hamiltonian. While 
these methods correctly include many-body effects they 
often require some knowledge of the empty state density 
of states (DOS) which cannot be easily obtained without 
the use of ab initio electronic structure calculations. For 
this reason core level spectroscopy is an ideal field for the 
application of combined ab initio and many-body cluster 
model approaches. 

The orthorhombic structure of La2Cu04 with the CuO 
planes oriented perpendicular to the c-axis allows the 
incident X-ray polarization to be unequivocally defined. 
Our twinned La2Cu04 single crystal was a plate-like rect- 
angular sample with the c-axis perpendicular to the sam- 
ple plane and the a/b axis parallel to the sides of the sam- 
ple. The experiment was performed at beamlinc ID16 at 
the European Synchrotron Radiation Facility. The ex- 
perimental setup was that of a Rowland circle spectrom- 
eter with a horizontal scattering plane, using a spheri- 
cally focusing, bent Si(444) analyzer £|. The absorption 
spectrum of La2Cu04 was obtained by measuring the 
total fluorescence yield. We also measured the partial 
fluorescence yield spectrum by using the spectrometer 
set to the peak of the Ka± fluorescence line and varying 
the incident energy as for an absorption measurement. 
These were measured for two orientations of the sam- 
ple: In the first geometry (noted a/b henceforth) the a-b 
plane was in the horizontal scattering plane such that it 
contained both the incident photon momentum and the 
polarization. They were both oriented at 45° with re- 
spect to the a/b axis. In the second geometry (noted c) 
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FIG. 1: Partial fluorescence yield absorption spectra for 
La2Cu04 in the a-b geometry (dashed line) showing tran- 
sitions to Cu 4p CT orbitals and the c geometry (solid line) 
showing transitions to Cu 4p,r orbitals. Each spectrum is 
normalized to its peak. The inset is a zoom of the pre-edge 
area where resonant scans are taken. 



the a-b plane, oriented vertically, was perpendicular to 
the scattering plane. The polarization was at 75° with 
respect to the a-b plane (15° with respect to the c-axis 
of the sample). These geometries were also used for the 
measurement of the ls2p resonant emission. It should be 
noted that due to the in-plane polarization of the photons 
the a/b geometry is sensitive to the in-plane or a orbitals 
only while the c geometry is sensitive primarily to the it 
orbitals. The absorption spectrum is characterized by a 
small pre-peak at the onset, usually attributed at least 
partly to on-site Cu ls-3d quadrupolar transitions. In 
Fig. ^ we show the partial fluorescence yield spectra for 
both geometries. In our RIXS measurements the inci- 
dent energy is limited to the region around the pre-peak 
(inset of Fig. [T]and left panels of Fig. Inland Fig. 0). It 
will be seen that the dipolar contribution in this pre-edge 
region being small, the weak on-site Cu ls-3d quadrupo- 
lar transition is easily visible. Most importantly we will 
show that an even weaker off-site Cu ls-3d dipolar tran- 
sition mediated by 02p states and never before identified 
is also clearly visible in the RIXS spectra. Though the 
absorption spectrum is relatively featureless, the emis- 
sion spectra are rich in structure and strongly dependent 
on the incident energy 

In Fig. [3and Fig. [3] (center panels) we show the experi- 
mental RIXS spectra for incident energies increasing from 
bottom up. The baseline of each RIXS spectrum corre- 
sponds to the incident energy at which it was measured 
as seen on the absorption spectrum on the left panel. 
RIXS spectra have been arbitrarily shifted such that the 
lowest energy transfer feature corresponds to the zero of 
the energy scale and have been normalized to the most 
intense feature in each spectrum. They show twin fea- 
tures separated by 20 eV due to spin-orbit splitting of the 
2p core level. In the following we shall discuss only the 
region from eV to -20eV energy transfer. The spectra 
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FIG. 2: a/b geometry. Left panel: The pre-edge region of 
the absorption spectrum. Experiment (dotted line), cluster- 
model calculated on-site Cu ls-3d quadrupolar (dashed line) 
and off-site Cu ls-3d dipolar contributions (solid line) . Center 
panel: Measured RIXS spectra. The incident energy can be 
traced to the absorption curve in the left panel. Right panel: 
Calculated RIXS spectra using ab initio Cu 4p DOS (solid 
line) and the modified DOS shown in the inset of Fig. [I] 
(dashed line). The peak at eV is due to the on-site Cu ls- 
3d transition. The additional peak appearing at -3 eV in the 
RIXS spectra is due to the lowest energy feature of the 4p CT 
DOS which is determined by the hybridization with off-site 
Cu 3d empty states. 




Transferred Energy (eV) Transferred Energy (eV) 

FIG. 3: c geometry. Left panel: The pre-edge region of 
the absorption spectrum. Experiment (dotted line), calcu- 
lated on-site Cu ls-3d quadrupolar contribution (dashed line). 
Center panel: Measured RIXS spectra. The incident energy 
can be traced to the absorption curve in the left panel. Right 
panel: Calculated RIXS spectra. 

are dominated by a sharp peak at eV which resonates 
when excitation energies are on the absorption pre-peak. 
In CuO a similar structure has been interpreted as due 
to a quadrupolar transition^, By separately cal- 
culating the contribution of Is — > 3d quadrupolar, and 
Is — > 4p a and Is — > Ap^ dipolar transitions we easily con- 
firm that in La2Cu04 this peak is due to the quadrupolar 
transition. In both geometries two additional features are 
seen. In the a/b geometry (Fig. |5J a weak, broad, peak 
is observed around -17 cV while a sharp peak at about 
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-3 eV resonates only at higher incident energies. In the 
c geometry (Fig. |3J) a strong broad structure with twin 
peaks at -12 eV and -7 eV is seen while a second peak 
appears at -5 eV, resonating at higher incident energies. 
The calculated results of the Cu Ka RIXS are shown in 
Fig. Inland Fig. (right panels) and are in excellent agree- 
ment with the experimental ones. The broad structures 
at higher energy loss are attributed to virtual transitions 
due to the main peaks in the Cu 4p DOS and are slightly 
shifted in energy from the measured ones. The lower en- 
ergy resonant structures are due to transitions to features 
in the tails of the DOS. These attributions were made by 
modifying portions of the calculated DOS and tracing 
the modifications in the calculated RIXS spectrum. 

RIXS data are interpreted using a unified ab initio 
and multiplet calculations approach. Numerical multi- 
plet calculations were carried out based on CuOg°~ sin- 
gle cluster model with intra-atomic full multiplets in D^ 
local symmetry and including the screening channel for 
charge transfer from O 2p ligand state. Only a single Cu 
atom is retained and charge transfer between the Cu 3d 
state and the neighboring O 2p orbitals is allowed. The 
Hamiltonian is given by H = H + H pd + H ch + H nm i tip i ct . 
H includes the on-site energies of the Cu Is (£i s ), 
Cu 2p (e 2p ), Cu 3d (e 3d (T),), Cu 4p (e 4p (k)) and Lig- 
and O 2p (e p (r)) states; H = £ r>(r £3d(r)4vdiv + 

J2r a £ p(T) a Ta a ro-- The irreducible representations of 
the D 4 /j symmetry are T ( = b\ g , a lg , b 2g and e g ). H pd 
is the hybridization between Cu 3d and O 2p ligand 
states, H pd = J2r,v ^( r )(<4v a r<r + al a d ra ). The core- 
hole interaction H dc is composed of two contributions, 
the attraction between Cu 3d electrons and Cu Is holes 
(Ud c (ls)) an d between Cu 3d electrons and Cu 2p holes, 
(Udc(2p)), namely H dc = U dc (ls) Y,v,<r,a> 4v dr v(l - 

sls a/ )+U dc (2p)J2r., m ,a,a' 4<A<7 (l-p^JW ) • Finally 
the Hamiltonian i/ mu itipiot describes the intra-atomic 
multiplet coupling originating from the multipolc compo- 
nents of the Coulomb interaction between Cu 3d states 
and that between Cu 3d and 2p or Is states and the spin- 
orbit interactions for Cu 3d and 2p states 0, 0, 0] • 

The parameter values used are (in eV): Ud c (ls) — 
Udc(2p) = -7.0, V(b lg ) = 2.5, V(a lg ) = 1.44, V(b 2g ) = 
1.25, V(e g ) = 0.88, A = 2.75(charge transfer energy), 
T pp = 1.0 (hybridization between nearest neighbor O 2p 
orbitals), very close to the ones used in ref. 6]. The 
Slater integrals and the spin-orbit coupling constants are 
calculated by Cowan's Hartree-Fock program, 0] and 
then the Slater integrals are scaled down to 80%. 

The Ka RIXS cross-section is calculated on the basis of 
the formula of the coherent second order optical process 
as S(ft,oj) — S£>(Q,ij) + Sq(CI,u>) where S a (Cl,u>) with 



a = D,Q is: 
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and g), \ m) and | /) are the ground, intermediate and 
final states of the Hamiltonian H with energies E g ,E m 
and Ef, respectively. The incident and emitted photon 
energies are represented by and w, respectively. The 
core-hole lifetime broadening is denoted by T K = 0.4 eV, 
half width at half maximum (HWHM) for the Is core- 
hole and T L = 0.2 eV (HWHM) for the 2p core-hole. 
Gaussian broadening due to the experimental resolution 
is taken to be 0.6 eV (HWHM). The operators T a with 
a = D,Q represent the optical dipole and quadrupole 
transition respectively. The polarization of the incident 
photon is also taken into account III • 

To calculate eq. ^ we use > as basis states, two con- 
figurations, 3d 9 and 3d 10 L, where L represents a hole in 
the ligand states. The basis configurations in the inter- 
mediate state | m) are ls3d 10 for the quadrupole tran- 
sition and ls3d 9 4jo and ls3d 10 IAp for the dipole excita- 
tion. The final states | /) are 2p3d 10 (quadrupole) and 
2p3d 9 4p and 2p3d 10 L4p (dipole). In the case of dipolar 
transitions we assume the intermediate states | m) (final 
states | /)) to be a direct product of a Cu 4p state | k) 
and the remaining system | m') ( | /') =| k) \ m')) with 
energies E m = E m > + Eu (Ef = Ef + £fc). Thus we can 
transform the sum over m in eq. \I]to a sum over m' plus 
an integral over the high enenergy continuum weighted 
over the 4p Cu projected density of states p v {e) (with 
77 = cr,7r), namely: J2 m i — ► J2 m > I de Pv( e m- 

The Cu 4p projected density of states, p v (e) = 
j^Z kX fen l< - 0|Cu i 4ff T? )| 2 (5(e kn - e), is calculated ab 
initio^ using the Kohn-Sham eigenvalues (ekn) and 
eigenfunctions (ipkn) and the |Cu'4p^) orthonormalized 
atomic orbitals. We disregard the core hole potential act- 
ing on the 4p state, because it is screened strongly by a 3d 
electron transferred from the O 2p states. Since we con- 
sider a supercell with two La 2 Cu04 molecules, p v is also 
summed over the two inequivalent Cu atoms (i = 1,2). 
In the ab initio simulations we use ultrasoft pseudopo- 
tentials 01 w ith a 45 Ry cutoff and the spin-polarized 
generalized gradient approximation(SPGGA) [llj for the 
exchange-correlation kernel. Since using the SPGGA 
the ground state turns out to be metallic we adopted 
the SPGGA+U approximation [T3| with U = 8eV in 
the implementation of ref. Self-consistence is per- 

formed over a 6 x 6 x 6 k-points mesh in the Brillouin 
zone centered at the T point. We converge to an an- 
tiferromagnetic state with Cu local magnetic moment 
0.625pB, to be compared with the experimental value 
of 0.55 ~ O.6O/is0,^|- The band gap is approximately 
0.65 eV, in good agreement with previous calculations of 
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FIG. 4: Projected Cu 4p^ (solid) and 4p a (dashed line) or- 
bital density of states calculated ab initio . The inset shows 
the low energy region in detail. The hatched peak, the lowest 
energy structure in the Cu 4p a DOS, is induced by hydridiza- 
tion between the Cu 4pa:, H atomic orbitals, the empty O 2p 
orbitals and the empty 3d orbitals of the neighbouring Cu 
atoms forming a 6i s state. 



ref. (0.7 eV), but substantially smaller than the one 
determined from optical conductivity data in ref. 01 
(rs 2 eV). The Cu 4p CTi7r projected DOS calculated on a 
Nk — 20 x 20 x 20 k-points mesh are shown in Fig. ^ 

For a finer comparison between theory and experiment 
let us first consider the sharp peak at eV in Fig. [21 
and Fig. due to the on-site ls-3d quadrupolar excita- 
tions, where the intermediate state is ls3d 10 and the final 
state is 2p3<i 10 . The 2p3d 10 final state shows no multiplct 
structure because the 3d shell is filled, which explains the 
single peak spectral shape. Next, we consider the dipolar 
contribution. In the a/b geometry the polarization being 
in the a-b plane, only Is — > Ap a transitions contribute. In 
the c geometry n transitions are visible, the a component 
being almost negligible (see Fig. |3J) • Finally we illustrate 
the extreme sensitivity of this method by looking for the 
origin of the peak at -3 eV in the a/b geometry. A closer 
look at the Cu 4p DOS in Fig. 0]reveals that though the 
main peak of the a component is at higher energy than 
that of the n component, the lowest energy structure is a 
in nature. The intensity of this structure, peaking at leV, 
is less than one percent of the total DOS in this region. 
Nevertheless if this structure is artificially truncated (in- 
dicated by the hatching in Fig. g), the peak at -3 eV 
disappears in the calculated RIXS (Fig. [21 right panel, 
dashed curve). This low energy structure in the Cu 4p CT 
DOS is induced by the hybridization with O 2p states 
and with the Cu 3d b\ g empty state of the neighbouring 
Cu sites. This is the first identification of this ls-3d off- 
site transition mediated by the O 2p orbitals and shows 
the chemical detail which such a measurement reveals. It 
is also present in the calculated absorption in the a/b ge- 
ometry (Fig. left panel) while absent in the calculated 
absorption in the c geometry (Fig. [21 left panel) and in 
the measured and calculated RIXS (Fig. [31 center and 



right panels) 

XANES measurements have long been the standard 
method for accessing excited states. Earlier works 
[lH lifT ] have pointed to the advantages of RIXS. We 
have conclusively shown with the example of La2Cu04 
that RIXS at the Cu Ka line provides a sensitive and 
precise way to determine the Cu 4p unoccupied DOS 
and the interactions of the Cu atom with its immediate 
environment. We distinguish between the dipolar and 
quadrupolar part of the transitions and separate them 
according to their symmetry. We also find that a new 
approach combining ab initio calculations of the partial 
Cu-4p DOS with many body cluster calculations is very 
successful in explaining the data. Indeed both the po- 
larization and incident energy dependance of the spectra 
are fully explained. Features in the measured spectra are 
traced to fine details of the calculated Cu-4p DOS (used 
directly) giving a very sensitive method for checking its 
reliability and identifying for the first time the off-site 
Cu ls-3d dipolar transition. In particular this work, by 
validating the nature of the Cu 4p unoccupied DOS, will 
lead to the definitive attribution of the origin of the differ- 
ent features in the Cu K-edge X-ray absorption spectrum 
[2ll |22| . Finally this method should prove useful for in- 
vestigating the change in these excited states produced 
by doping as well as those induced by changes in local 
symmetry, notably those between the Lanthanum based 
hole-doped compounds and Neodymium based electron 
doped compounds. M.C. acknowledges discussions with 
G. A. Sawatzky. This work benefitted from CNRS/JSPS 
research grant No. 16030. Calculations were performed 
at the IDRIS super-computing center (grant no. 051202). 
SWC was supported by the NSF-DMR-0405682. 
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